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NOMENCLATURE 

K housing scale factor 

KR housing radius 

L seal length 

p pressure 

R shaft radius 

r radius 

T torque 

V velocity 

x void fraction or void region 

x,i local coordinate 

Y parameter; see eq. (20) 

Z parameter; see eq. (12) 

a operator, see eqs. (8) to (11) 

u viscosity 

p density 

t shear 

to angular velocity 

Subscripts: 

g gas 

£ liquid 

m limiting rotation boundary 

o pseudo boundary, equivalent to KR 
r radial 

ref reference boundary 

e angular 

INTRODUCTION 

Two-phase flow in shaft seals, bearings, pumps, 
and dampers has been of concern to tribologists for a 
number of years (1). The void phase can be either 
vapor or dissolved gases (inerts) and in either case, 
a bubbly lubricating film can occur that may substan- 
tially affect performance. 

Several investigators have addressed the favorable 
effects of bubbly films in steadily loaded seals and 
bearings (e.g., 2-4). However, under dynamical loading 
with squeeze effects, there appears to be a substantial 
reduction in the load capacity of bearings (5,6); the 


effect in seals is still under investigation. In most 
instances one trades load capacity (stiffness) for 
damping. 

The entrainment of void (air) into oil is common- 
place and significant cavity regions can form in jour- 
nal bearings (7). It is clear that body force effects 
due to rotation or cyclic motions along with continuity 
requirements have a significant effect on void migra- 
tion and coalescence. 

In this study air-entrained oil was discovered to 
migrate toward and in some instances become attached 
to the rotating interface. A Bingham fluid was used to 
model the phenomena. 

APPARATUS AND INSTRUMENTATION 

The apparatus consists of a 16-mm (5/8-in.) diam- 
eter shaft rotating at 1650 rpm, approximately 82.5 mm 
(3-1/4 in.) from the bottom of a 230- by 230- by 350-mm 
(9- by 9- by 14-in.) oil reservoir (Fig. 1). A flex- 
ible, transparent cylinder could also be used to de- 
crease the effective size of the reservoir to that 
approaching a seal or bearing. Rotational speed varied 
in a continuous manner. In some tests the shaft con- 
tainer interface was sealed by cooling the oil to form 
a natural plug (8). The entrained bubbles were formed 
simply by pouring oil into the reservoir. The en- 
trained bubble motion and shaft rotation were photo- 
graphed at 400 pictures/sec. 

ANALYSIS 

It is well known from balancing the radial forces 
that an entrained gas in a liquid will migrate against 
the body force field and establish itself at a balance 
point - which, in this case, is either about the cen- 
terline or about the shaft. The puzzling part is that 
the bubble appears as nearly solid-body rotation - the 
bubbles approach one shaft radius in diameter and are 
whirling about with the oil. However, if one estab- 
lishes the flow field for a shaft (r = R) in an infi- 
nite reservoir, the flow field is nearly solid body 
close to the shaft (R < r < KR). Apparently, the 



outer limit of the oil flow field possesses a shear 
layer that, not only limits the escape of the vapor, 
but distorts the bubble as well (Fig. 2). 



Single Phase 

For a single-phase fluid bounded between two con- 
centric cylinders, the torque applied by the rotating 
inner cylinder (shaft) on the fluid is found from a 
solution to the momentum equation 

T -- <1^4 (1 
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where 


t:t < 2t 
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and the rate of deformation tensor is 


where 


<T) = 4.UR 2 



Here, R is the inner cylinder radius, KR is the outer 
cylinder radius, and the angular velocity is 


= t:aV - | (v*V ) 


( 8 ) 


(9) 


( 10 ) 



r K 2 R 
R " r 


which is commonly known as the Bingham plastic model. 
For our model t re is the dominant term 


It is clear that the ratio of the torque applied to a 
gas to that applied to a liquid is determined by the 
ratio of viscosities 
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and solving gives the circumferential velocity as 


fl 


R < r < r 
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Two-Phase 

If it were possible to segregate the fluid domain 
encompassing the shaft into x domains of vapor (gas) 
and (1 - x) domains of liquid, assuming a 1 iquid/ vapor- 
slip ratio of unity, the torque would become 




where r Q = KR, 

Z - •jr- (13) 
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However, eq. (5) possesses a singularity at the satura- 
ted vapor locus and a better representation becomes 
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where 


<h> ■ <0 


(7) 


From eq. (6), it is clear that the torque is signifi- 
cantly reduced by adding a small amount of gas (vapor). 


and 


< T o) = 


T 


(14) 


Further, assuming that KR = r 0 is known, equa- 
ting the torque applied to the vapor between the con- 
centric cylinders to that of the Bingham model provides 
an expression for t 0 . 



(15) 


where 


Attached Bubble 


For a bubble attached to the rotating shaft, we 
assume that the flow field and outer cylinder can be 
replaced by a fluid with a stress tensor defined by 
eq. (8). The outer cylinder boundary KR is now re- 
placed by a pseudo boundary r 0 . 


K 


o 


R 


(16) 


The circumferential velocity of eqs. (12) and (13) be- 
comes 


2 




where r m is the radius when V e = 0. 

To establish a limiting value for r m , V e is set 
equal to 0, and eq. (17) becomes 



and 

= 2 - Y' 1 + In Y (19) 

where 


tures taken at 400 frames/sec revealed individual gas 
cavities adhering to the shaft and rotating at nearly 
synchronous speed (Fig. 3). The bubbles were rather 
distorted at their apex because of shear at the inter- 
face and appeared to locally "lock“ oil between them 
to complete the "solid-body rotation." These bubbles 
extended 0.6 shaft radii from the shaft surface; they 
also tended to migrate along the shaft and to coa- 
lesce. Large bubbles would lose gas to the oil. That 
gas would then appear as tracers in the oil. When the 
seals were not cooled, some void would migrate toward 
the seal aperture and disrupt the flow (reduce the 
leakage) . 

Applications to other data for seals, bearings, 
and pumps have not been addressed, but such entrainment 
effects should be expected when the fluid is two-phase. 
At higher rotational speeds and during transition to 
turbulence, the bubble size and rotational synchroniza- 
tion are expected to be limited. Further, as the res- 
ervoir approaches that for the classic flows in rotat- 
ing annuli, the entrained bubbles are expected to merge 
into Taylor vortices. 


Y = 



( 20 ) 


For attached bubbles with diameters of 0.5, 1, and 2 
shaft radii, the values of r m are approximately 25 mm 
(2 in.), 130 mm (50 in.), and • . 

With such a definition of the shear layer, one can 
expect nearly solid-body rotation to a radius of r 0 
with viscous flow at larger values of r. Since the 
radial pressure gradient is given directly in terms of 
the angular velocity. 


dp 

Hr 


oV‘ 


(21) 


it is clear that the vapor is 'held' adjacent to the 
shaft by the difference in densities, or 


dp 
Hr : 




( 22 ) 


and that this is strongly dependent on the rotational 
speed. Further, as the interface r = r 0 is ap- 
proached, the gas entrained into a bubble should become 
distorted by the change in shear as defined by eq. (12) 
or (17). 

EXPERIMENTAL RESULTS 

To achieve gas or vapor entrainment into the 
fluid, an oil was poured into the reservoir, and the 
shaft rotated at selected speed. The vapor adjacent 
to the shaft rapidly formed what appeared to the un- 
aided eye as a vapor torus with a slower fluid/bubble 
motion within the reservoir. High-speed motion pic- 


SUMMARY 

An experimental program has been carried out to 
determine the nature of the flow field about entrained 
gases in liquids (air/oil in this study). The en- 
trained gas formed bubbles that were observed to rotate 
with the shaft in nearly a synchronous manner for se- 
lected rotational speeds. A simple analysis using a 
Bingham fluid was used to describe the phenomena. 
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Fig. 1. - Schematic of experimental apparatus. 
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Fig. 2. - Model of bubble entrainment about a 
rotating shaft. 



(a) PLAN VIEW (b) EDGE VIEW 


Fig. 3. - Entrainment of air bubbles in oil about a rotating shaft. 
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